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ABSTRACT 


The Doppler factors occurring in range- rate satellite-to- 
satellite tracking measurements are derived with special 
relativistic effects included. The error resulting from 
the use of simplified expressions for these factors is 
discussed. 
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DOPPLER FACTORS IN SATE LLITE-TO- SATELLITE TRACKING 


INTRODUCTION 

The Doppler factors arising in one-way and two-way satellite tracking are well- 
known [1, 2] . The factors applicable to the satellite- to-satellite tracking con- 
figuration are not so readily available, however, and a special relativistic 
derivation of these factors is documented in this report. 


THE DOPPLER EFFECT FORMULA 


The relativistic formula for the Doppler effect [3] is often given in the form 


y* 1 - v 2 /c 2 
1 - e ♦ v/c 


(1) 


In (1) f Q is the proper frequency of the source (i.e. the frequency measured in 
the inertial system in which the source is at rest). The source itself is moving 
at a vector velocity v with respect to an observer. The frequency detected by a 
receiver at rest in the inertial system of the observer is f. The vector e is a 
unit vector in the direction of propagation as measured in the coordinate system 
of the observer, v is the magnitude of the vector v, and c is the speed of light. 

If the distance between the vehicle and the observer at the time of emission of 
the signal is r, then (letting the vector position of the vehicle be f , differenti- 
ating ? • r = r 2 with respect to time, and setting ? = -re) 


e 


v - 



and Equation (1) provides the downlink Doppler factor 


f/f 0 


l/l - vVc’ 
1 + r/c 


( 2 ) 


( 3 ) 


that holds good for one-way range-rate tracking, the Doppler factor being here 


1 



defined as the factor by which the transmitted frequency is multiplied to give 
the received frequency. To obtain the Doppler factors applicable in satellite-to- 
sateliite tracking, however, a more general form of the Doppler formula (1) is 
required; one in which not only the transmitter but also the receiver of the signal 
is in motion with respect to the observer. 


GENERAL DOPPLER EFFECT FORMULA 

The required generalization of the Doppler formula (1) can be obtained from the 
latter by use of standard transformations of special relativity. Consider an ob- 
server at rest in an inertial frame I, a transmitter moving with vector velocity 
u in S, and a receiver imbedded in a frame 2' that moves with vector velocity v 
with respect to the observer. The configuration is depicted in Figure 1. As 



Figure 1. Doppler Effect Geometry 
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seen by the observer in 2, let i be the displacement vector drawn from the 
point a where the signal is emitted at time t a to the point b where the signal is 
received at time t . Then the unit vector e in the direction of propagation is 


e/e 


( 4 ) 


where e is the magnitude of e and is equal to 


e = c(t b - t a ) (5) 

Following convention, the quantities above, being measured by the observer in 2, 
are unprimed. The same quantities, as measured by an observer in 2', will be 
identified by a prime. With this notation, the Doppler Equation (1), as it applies 
to an observer in 2' becomes 


f' 


V 1 - u' 2 /c 2 
1 - e' • u '/c 


( 6 ) 


Here u' is the magnitude in 2' of the vector velocity u' of the transmitter. 

To obtain the desired form of the Doppler equation, the quantities u', u\ and e' 
in (6) must be expressed in terms of the unprimed quantities seen by the observer 
in 2. The vector equation for the relativistic transformation of velocities is 
[3, Eqn. 2.55] 


yi - v 2 /c 2 u + [( i _ yr ~ v 2 /c 2 ) u ■ v/v 2 - ij v 


1 “ u • v/c 2 


(?) 


Taking the scalar product of (7) with itself, and performing some algebraic 
manipulations, the well-known result [3, Eqn. 2.56] 



yi - v 2 /c 2 y i - u 2 /c 2 

1 - u • v/ c 2 


( 8 ) 


is obtained. 
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To express e' in terms of unprimed quantities, the Lorentz transformation in 
vector form [3, Eqn. 2.25] is used. 


x -I- V 



( i / yi - v 2 /c 2 


i) - 1/ yi - v 2 /c 2 


(9a) 



(9b) 


Here (x, t) are the coordinates of any given event. Substituting (b, t b ) and (a, t a ) 
successively into (9a), taking the difference of the two resulting equations, and 
employing (5) 


— t 

e 


e + v 


i — 

e ‘ 

v : 


(1/ y 1 - v 2 / c 2 - 1) - e/(c "j/i - v 2 7c 2 ) 


( 10 ) 


Substituting in the same way into (9b), and multiplying the resulting difference 
by c, 


e' - (e ~ v • e/c)/ yi - v 2 /c 2 (11) 

Taking the scalar product of (7) and (10), dividing both by (11) and by the constant 
c, and subtracting the result from unity, there obtains 

- (/ _ (1 - e • u/c) (1 ' v 2 /c 2 ) 

1 - e • u /c - — — ^ 

(l~e • v/c) (l - u ' v/c 2 ) 

Substituting (8) and (12) into (6) 

(1 - e ■ v/c) V 1 - uVc 2 
f' = f 0 ~-r 

(1 - e * u/c) yi - v 2 /c 2 


( 12 ) 


(13) 
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Equation (13) is the desired generalized Doppler formula. It relates the fre- 
quency f 0 measured aboard a transmitting vehicle moving with vector velocity u 
at the time of transmission of a signal, to the frequency received aboard a 
vehicle moving with vector velocity v, at the time of reception of the signal the 
unit vector e being directed from the point where the transmitter is located at 
the time of transmission to the point where the receiver is located at the time 
of reception. Using this equation all of the Doppler factors involved in satellite- 
to-satellite tracking can be written down immediately by suitably specifying the 
transmitter and receiver velocities. 

Note, incidentally, that (13) may be written in the form 


(f']/l - vVc 2 ) 


(f 0 yr^W) 


1 - e • v/c 
1 - e • u/c 


(14) 


from which it becomes apparent that the square roots that appear in (13) may be 
regarded as corrections to the classical Doppler factor (1 - e ♦ v/c) / (1 - e • u/c) 
[3, Eqn. 1.22] for the effect of relativistic time dilation on the frequencies mea- 
sured aboard the vehicles. 


SATELLITE-TO-SATELLITE TRACKING DOPPLER FACTORS 


Referring to Figure 2, let a signal from a ground-based tracking station be re- 
ceived at time t l by a first satellite moving with a velocity v. Let s be the vector 
drawn from the tracking station to the point where the signal is received. After 
an interval corresponding to the transponder delay, a signal is sent from the 
first satellite at time tj to a second satellite where it is received at time t 3 . 

Let q be the displacement vector drawn from the point of transmission at time t 2 
to the point of reception at time t 3 . In a similar manner, let the round-trip of 
the signal be completed by retransmission back to the ground station along the 
vector paths p and r. Note, however, that the direction of the latter vectors has 
been taken opposite to the direction of propagation. 

To obtain the uplink Doppler factor u,, the quantities u and v appearing in (13) 
are defined to be the vector velocities of the tracking station and the first satel- 
lite respectively. Thus in (13) u = 0 and e = s/s = s giving 

1 - s /c 

“i " ■ ; (15) 

yi - v 2 (s) / c 2 
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where v(s) is the speed of the first satellite at the time t x when the signal was 
received at the vector position s. 

The downlink Doppler factor d x has already been given in (3) 


l/l - v 2 ( r ) / c 2 

d x = — 7 

1 + r/c 


(16) 


The second uplink Doppler factor u 2 is obtained from (13) by interchanging u 
and v 


u 


2 


(1 - Ug / c ) }/l - V 2 (q)/ c 2 

(1 ' V - /c) yi - u =(q)7? 


(17) 


Here u q = q • u is the component of the velocity of the second satellite along the 
q direction, and v q is the like component of the first satellite. 

The remaining Doppler factor d 2 follows from (13) by setting e = - p giving 


d 2 


(1 + v p / c ) |/i - u 2 (p)/c 2 
(1 + u p /c) yi - v 2 (p)/c 2 


(18) 


In coherent round-trip Doppler tracking, the Doppler factors appear only as 
products, for example u x d x or u x u 2 d x d 2 . In such cases, the square roots ap- 
pearing in (15)- (18) cancel out provided that the speed of the vehicles does not 
change appreciably over the duration of the round-trip. Consequently negligible 
error is incurred in these cases if the classical Doppler factors 


1 - s/c 

(19) 

1/(1 + r/c) 

( 20 ) 


,7 



( 21 ) 


U 


2 


1 ~ U q/ C 

1 “ v q /c 


1 + v p /c 
d 2 ' 1 + U p / C 


( 22 ) 


are used in place of (15)-(18). 


APPROXIMATE DOPPLER FACTORS 

If the first satellite in Figure 2 is geo stationery, which is a case of practical 
interest [4] , it is possible to further simplify (21) and (22) because of the rela- 
tively small magnitude of the velocity of the geostationary satellite. The de- 
nominator of (21) can be expanded to give 


u 


2 


(1 - u q /c) 





= 1 - 7 (1 + v q /c) 

since u q - v q = q. Consequently the approximation 

u 2 - 1 - q/c 

produces a relative error e in range rate that is less than 

V q 

e T 1 — 1 V/C 


(23) 


(24) 


( 25 ) 
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and an absolute error less than 


e a < qv/c 


(26) 


Typical values for u and v in the ATS/NIMBUS-F experiment are ( 5] 7.5 km/sec 
and 0.08 km/sec giving 


e a < 2.0 millimeters/sec (27) 

An error of approximately the same magnitude but of opposite sign arises in 
using 


d 


2 


1 


1 + p/c 


(28) 


in place of (22) . 

An error of 2.0 mm/ sec is not insignificant. However this error, being of oppo- 
site sign in (24) and (28), cancels out in round-trip Doppler measurements, in 
which only the products of uplink and downlink Doppler factors appear. A more 
complete discussion of the errors is given in the Appendix. 


SUMMARY 

Using the transformations of special relativity, the four Doppler factors (15), 
(16), (17), and (18) encountered in satellite-to-satellite tracking are derived. It 
is shown, in the case of round-trip tracking from the ground through a geosta- 
tionary satellite to a close-in earth orbiting satellite and then back to the ground 
view the same path, that these Doppler factors may be replaced by the simpler 
ones given by (19), (20) and either (21) and (22) or (24) and (28). 
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APPENDIX 


APPROXIMATION ERRORS 
A. NEGLECT OF RELATIVISTIC FACTORS 

Using Equations (15) and (16) the quantity c (1 - Ujd^/2 may be expanded in in- 
verse powers of c to give 


c (1 - Uj dj) / 2 


r + s r (r + s) v 2 (r) - v 2 (s) 
2 2c + 4c 


(Al) 


where terms in negative powers of c other than the first power have been neg- 
lected. If Equations (19) and (20) are used instead of (15) and (16) to calculate 
c (1 - u 1 d 1 )/2 there results ' 


c (1 


u i-di)/2 



r (r t s) 


2 c 


+ 


(A2) 


The first term on the right hand side of (Al) is the mean range-rate of the first 
satellite, and this term is much larger than the remaining ones. It follows that 
the use of the classical Equations (19) and (20) instead of relativistic (15) and (16) 
will result in an error in the calculation of range-rate of the magnitude 


v 2 (r) ~ v 2 (s) 
4 c 


(A3) 


The numerator of (A3) represents the change in the square of the speed of the 
satellite during the time it travels from point s to r. 


1 

4 c 


Av 2 = 


1 _ 
2 c 


d v 

v dt At 


(A4) 


In the case of an earth- orbiting satellite in unpowered flight, dv/dt will 
not exceed the acceleration of gravity at the earth's surface (about 9.8 m/s 2 ). 
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In single satellite tracking the time interval in question is a transponder delay 
that might be, very roughly, about 1 microsecond. Using these values, and 
7.5 km/s for v there results 


6 a _< 1. 2 x 10 -10 m/sec 


(A5) 


which is clearly negligible. Like results will obtain in satellite-to-satellite 
tracking where the cancellation of relativistic factors in the product u t u 2 di d 2 
takes place between factors separated by successive transponder delays (for 
example, the radical in u 2 is cancelled not by the radical in d t but rather by one 
of the radicals in u 2 ). 


B. USE OF SIMPLIFIED FACTORS 

Expansion of the c (1 - u 2 d 2 )/2 using (21) and (22) gives 

c (1 - u 2 d 2 )/2 = -y (q + p) ^1 - ~j - 2 “ (v p p - Vq q) ( A6 ) 

Using the approximate Doppler factors (24) and (28), instead of (21) and (22) 
yields 


c(l-u 2 d 2 )/2 - -^(q + p) ( A7 > 

The error incurred in using (24) and (28) to calculate the mean range rate (q + p)/2 
between the two satellites is therefore 


€ a “ OpP “ v q q)/2c 


(A8) 


which is half the difference between the errors (see 26) incurred in using (28) 
or (24) separately. The error here is the increment 
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where Ar = t 4 - t 3 is the transponder delay in the second satellite and At = 
t s - t 2 is given by 


At 


q + P 
c 


+ A r 


2 q/ c 


(A13) 


The size of the change in velocity of the second satellite is less than 


| A u. | < g A r 


(A14) 


where g = 9.8 m/s. 

Substituting (A12) and (414) into (All) 

A u < g A r + — ~ At + A t (A15) 

q " 6 q q 

Evaluating Av q in the same way 

a 

A v < g A t + — A t + — — A t (A16) 

q — B q q 


where here g may be taken to have a smaller value because of the remote loca- 
tion of the geostationary satellite. From (A15), (A16), a bound on e a , Equa- 
tion (A9) may be determined. The contribution from the first term on the right 
hand side of (A16) is much larger than all the rest, giving the result 

€ a 5 ugAt/2c (A17) 


Using u - 7.5 km/s, g = 0.2 m/s 2 , and At = 0.25 s 


e < 6.2 x 10- 7 m/s 

a 


which is negligible. 
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